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We compute the electric conductivity of nanowires in the presence of magnetic domain walls 
by the method of wave packet propagation. We demonstrate that the propagation through the 
wire depends on the initial state used in the wave packet simulation. We propose a procedure, 
based on the Landauer formula, to reduce this dependence. Direct numerical calculations of 
the Kubo formula for small individual systems are used as reference data for the proposed 
procedure. The enhancement of the conductivity in the magnetic domain wall is found in the 
present method in accordance with the previous perturbative analysis. 

KEYWORDS: Electric conductivity in magnetic domain wall, time dependent Schrodinger equation, Kubo for- 
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§1. Introduction 

In general, transport phenomena are intimately related to the dynamics of particles involved. In 
the case where the direct interaction between the particles that participate in the transport process 
are negligible, the solution of of the one-particle time-dependent Schrodinger equation contains 
all information, at least in principle. Of course, the same information is obtained by solving the 
corresponding stationary problem. Appealing features of the time dependent approach are its 
conceptual simplicity, its flexibility in terms of system geometry and choice of (spin-dependent) 
scattering potentials, and its absolute numerical stability. On the other hand, in contrast to the 
time independent approach, it is not well established how to extract the transport coefficients from 
the time dependent data. 

Recently, the electric conductance in nanoscale wires has attracted much interest. 1 ) In particular, 
the effect of the magnetic domain wall in ferromagnetic substances on the conductance has been 
studied extensively. The effect of a magnetic domain wall on the conductivity was clearly observed 
in experiments on nanowires. 2 ~ 7 ) In these experiments, the field is initially applied along the wire 
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direction and the magnetization uniformly aligns ferromagnetically in this direction, where the 
conductance is that of the wire with uniform magnetization. If the magnetic field is reversed, then 
the anti-domain of the magnetization appears at the edge of the wire. The appearance of the 
anti-domain creates a domain wall. The applied field is not very strong and hence only one domain 
wall exists. The anti-domain grows in time and finally the total magnetization is reversed. The 
resistance decreases shortly after the field is reversed and after some period returns to its initial 
value. This period corresponds to the life time of the domain wall. 6 ' 7 ) 

This phenomenon has been analyzed theoretically by various methods. 8 ~ 16 - ) One approach is to 
use the Kubo formula. Tatara and Fukuyama 8 ) pointed out that the existence of the domain wall 
causes the decrease of the resistance due to the quantum effect, which is in contradiction with the 
intuitive, classical picture that the current would be scattered with the domain wall. The key point 
is that the system has impurities which provide the resistivity of the medium and the domain wall 
suppresses the scattering by the impurities when the spins of electrons change adiabatically with 
the change of the magnetization. 

It is interesting to observe how the electron propagate through the magnetic domain and how 
it is scattered microscopically. This physical picture was investigated by a direct method, i.e. 
by a direct simulation of electron wave packet propagation in a system with impurities and the 
magnetic domain wall. 9 ) We will call this approach 'wave packet method'. This direct method 
provides intuitively understandable information about the electron propagation. 

We could also extract information for the transport coefficient from the simulation. The conduc- 
tance is naively related to the amount of the probability to go through the domain, which can be 
regarded as the transmission coeffient. This quantity was studied to investigate the above problem 
of the conductance of nanowires with the domain wall, and the result supports the enhancement 
of the conductance (decrease of the resistance). 9 ) In this work, the conductivity is estimated by 
studying the propagation of a simple wave packet. However, the transmission coefficient is usually 
defined for the incoming plane wave instead of a wave packet. Actually, in the wave packet method, 
we have to specify the initial configuration of the wave packet. 

As shown below, we find that the transmission coefficient depends on the initial configuration 
of the packet. Therefore it is difficult to extract the conductivity from data obtained in a single 
initial configuration. Not only the quantitative amount but also qualitative features change with 
the initial configuration of the packet. Thus, we study how much information for the transport 
coefficient we can obtain by the wave packet method. 

In order to obtain the reference data for given systems we calculate the conductivity by direct 
numerical evaluation of the Kubo formula. In this paper, we mainly study the dependence of the 
conductivity on the width of the domain wall. We find that, in the case without the magnetic 
domain walls, the transmission coefficients obtained by the wave packet method with various ini- 
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tial configurations show common qualitative features. In this case we do not need an additional 
procedure to extract the qualitative features of the conductivity from the simulation data. 

However, as mentioned above, when the domain wall is present, the transmission coefficient 
strongly depends on the initial configuration, e.g., various shapes of wave packets and spin po- 
larizations. Under this circumstance, we need an additional procedure to obtain the conductance 
from the simulation data. Inspired by the Landauer formula, we propose a procedure to estimate 
the conductance from the transmission coefficient data for various initial configurations. 

We study in detail the dependence of the conductance on the width of the domain wall in some 
small systems for which the direct numerical calculation of the Kubo formula is feasible. In such 
small systems, different configurations of the random potential lead to large fluctuation of the 
dependence. Our procedure successfully reproduces the individual results. 

We also study the effect of the width of the wave packets in the direction of the propagation, 
discuss its relation to the damping factor in the Kubo formula and find that quantum interference 
plays an important role. 

Applying this method we also study the conductivity in the domain wall system. Qualitatively 
similar results are obtained by the wave packet method and the direct numerical evaluations of the 
Kubo formula. In this perturbative regime, our results agree with those of Tatara and Fukuyama. 8 ) 

§2. Time-Dependent Schrodinger Equation for Propagation of Wave Packet 

The electron propagation in the nanowire is described by the time-dependent Schrodinger equa- 
tion 17 ' 18 ) 

f) r fi 2 „ 

<M*)> (2-i) 



— V 2 + V(r) - y. B S ■ M(r) 



where cr(=t or j) is the spin of the electron, m is the mass, and V(r) denotes spin-independent 
impurity potential. Here M(r) is the magnetization of the medium, a is the Pauli matrix denoting 
the spin of the electron and a ■ M(r) represents the magnetic coupling between them. 

For numerical work, it is convenient to use dimensionless quantities. First we express the energy 
in unit of the Fermi energy Ep, a typical energy scale of the system. Using 

ft = r +' = 1 (o 2) 

h/^E^' h/Ep ' 1 ' } 

we have 

i^ = {-V' 2 + y'-a-M'}*'EHy, (2.3) 

where ^'(r^t') = ^(r, t), M' = M(r y ), V = V{r'). Hereafter, in our numerical work, we use 
renormalized quantities, i.e. we express the energy in units of Ep (e = E/E-p), distances in units 



of xq = H/\/2mEp = 1/kp = Xp/2ir (x' = x/xq), and momentum in units of po = hkp = \j2mEp 
ip' = p/po) an d drop the primes in the notation. Writing the wave function in the form of the 
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Eq.(2.1) reads 

For time evolution 



-V 2 + 



V \ / M z M x - iM y 



^(r,t + r) = e- iHT *l>(r,t ), 
we use the method of exponential decomposition up to the second order 



e x(A+B) = e xA/2 e xB e xA/2_ 



Here we take 



A 



2m V 



(2.4) 



(2.5) 



(2.6) 



(2.7) 



B = -n B d- M(r) + V{r), 
where at each position f, A and B are 2x2 matrices. For most calculations we used the second 
order decomposition and, as a check, we occasionaly used the fourth order algorithm. 22 ) 

In our numerical simulation, we use a real-space representation of the wave function. The wave 
function is given by its value at points of the lattice and the value of the spin. For example, in 
two dimensions, we divide the space into N x x N y cells and each cell is identified by the coordinate 
Thus the wave function can be represented as 

( ^t(M.t) \ 

V T (l,2,t) 



V'T (N x ,N y ,t) 

V>l(M,i) 
^l(l,2,t) 



(2- 



The wave function at cell (i, j) is given by ^ a (i,j,t). We use several formulae to approximate V 2 . 
For example, the standard 3-point formula in one direction reads 

d 2 f ^ fi-i,j - 2fi,j + fi+ij 



dx 2 



A 2 



+ 0(A d 



(2.9) 
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where A is a mesh size. We also use the 5-point formula 



d 2 f _ ~fi+2,j + 16/i + ij - 30/j,j + 16/j-ij - fj- 2 ,j 



dx 2 ~ ' ' 12A 2 



+ o(a 5 ), 



(2.10) 



and the 9-point formula 



d 2 f , 




> 



+ o(a 4 ). 



(2.11) 
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We used Eq.(2.11) in this paper. In order to avoid the reflection from the edges of the system, we 
used absorbing boundary conditions at both edges. As for the impurity potential, we randomly 
distribute the impurity sites where V(f) ^ with some concentration at each cell. 

The quantity of interest in this work is the conductance through the wire. In the wave packet 
method we compute the transmission coefficient, and convert it to the conductance from it. We 
determine the transmission coefficient by measuring the electron current through a virtual screen, 
as indicated in Fig.l. The total current in the x-direction is given by 



The transmission coefficient T is calculated from the amount of current through the right-most 
detection screen. Likewise the reflection coefficient R is calculated from the amount of current that 
goes through the left-most detection screen. 

In Fig. 2, we show an example of electron propagation. Initially, we prepare a wave packet moving 
from left to right with the group velocity. Then it is scattered by the impurities. Part of the wave is 
reflected and some part crosses the impurities region and is transmitted through the wire. Intensity 
arriving at the ends of the wire is being absorbed. Note that even after a fairly long time (see Fig. 2 
(f)), we find some intensity in the impurity region. 

§3. Relation between the transmission coefficient and the conductance 

3.1 Sensitive dependence on the initial configuration 

The conductance a is related to the transmission coefficient T through the Landauer formula 
a = T/(l — T). The transmission coefficient T is obtained by solving the scattering problem for 
an incoming plane wave. However, in the wave-packet method we cannot use the plane wave. 
Therefore, to obtain T, we solve the Schrodinger equation for a specific initial configuration, and 
study the dependence on the initial configurations, 

As initial configurations, we use the rippled-Gaussian wave-packet 




(2.12) 




(3.1) 
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with k y = 1,2,3 and 4, for both a =| and j. Thus we have 8 different initial configurations. The 
shapes of k y = 1,2,3 and 4 are depicted in Fig.3. For each n, we determine the value of k x such 



For the eight initial configurations, we compare the transmission coefficients in a system with 
one magnetic domain wall M(x,y) with some random positions of impurity V(x,y), 



where Mq, xq, and X w denote the height, the position of the center, and the width of the magnetic 
domain wall, respectively. If the transmission coefficient is insensitive to the initial configuration 
we may obtain information for the conductance from one simulation. In the presence of a mag- 
netic potential, however, we find that the transmission coefficient is very sensitive to the initial 
configuration. In Fig. 4 we show the dependence of the transmission coefficients on the width of a 
single domain wall. We find that each initial configuration results in very different transmission 
coefficients not only quantitatively but also qualitatively. That is, some increase with the width 
and others decrease. Thus we have to be careful to deduce properties of the conductance from the 
data of the transmission coefficients, and we need some statistical treatment. 

3.2 Comparison with the Kubo formula 

In order to check the results, we calculate the conductivity by the Kubo formula for the system 
with a given configuration of impurity potential and magnetization. The method is described in 
the Appendix. In this approach we need to obtain all eigenvalues and eigenvectors of the system, 
and therefore we can only treat small system up to 256 x 16. Although we can calculate much 
larger systems by the wave-packet method, in order to compare with results of the numerical study 
of the Kubo formula, we study the same 256 x 16 system by the wave-packet method. Here in the 
wave packet method, the region of 256 x 16 is assigned to the scattering region as shown in Fig. 1, 
and we attached the extra region of 208 x 16 as leads at both ends of scattering region. 

First we compare the results for the case without magnetic domain wall. In Fig.5, we show the 
dependence of the Drude formula (solid line) , and the formula (dashed line) that includes the weak 
localization correction obtained by a perturbative treatment of the Kubo formula 8 ) 



where n denotes the electron density, r is the scattering time (r _1 = 2imi\v \ 2 n/H, where rtj is the 
density of the impurities and \v\ is the strength of the impurity potential), I = {hk-pr/m) is the 
mean free path, and (L x , L y ) is size of the system. The first term in the right hand side of Eq.(3.4) 
corresponds to the Drude formula. In Fig.5 we plot the data obtained numerical study of the Kubo 



that (iP* g \H\Ag) = Ep. 




(3.2) 



(3.3) 




(3.4) 
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formula (closed circle) and those obtained by the wave packet method with n = 1, 2, 3 and 4 (cross, 
asterisk, open square, and closed square, respectively). These data do not depend on the spin. We 
find that all of them qualitatively agree with each other. 

3.3 Effects of magnetic potential 

When the magnetic potential is introduced, as we saw in Fig.4, the data vary strongly with the 
initial configurations. Thus we need some method to extract information of the conductance from 
the data. In order to obtain the conductivity by the Kubo formula, periodic boundary conditions are 
required. A single domain wall is not compatible with periodic boundary conditions. Therefore, for 
simplicity, we adopt a screw shaped magnetic potential to study the effect of magnetic scattering: 

tanh(^a) 

MJx, y) = M (3-5) 

V ,yJ cosh ^) 



[ tanh(- 



where M x 2 + M- 2 = Mq 2 , Mq, xq, and X w denote the height, the position of the center, and the 
width of the screw structure, respectively. The magnetization of the constant strength Mq rotates 
in the x — z plane. Here Mq is taken to be Q.2Ep. 

3.4 Transmission coefficient and conductivity 

Now we describe the statistical procedure that we use to process the numerical data obtained for 
different initial wave packets. In the spirit of the Landauer formula, different initial configurations 
may correspond to different channels. The transmission coefficients vary from channel to channel, 
and therefore the dependence of the transmission on the initial states appears in a natural way. 
Let Tij be the transmission coefficient from the mode i to j, and let Rij be the reflection coefficient 
from the mode i to j. Because we cannot distinguish between outgoing channels we sum up all the 
outgoing channels. We define the transmission and reflection coefficients for the incoming channel 
i by 

j j 
Following Azbel et al., 19 ) the conductance is given by 

^T^'^lS-T,)^' (3 - 8) 

where Vi = fiki/m. As to the incident modes, we take the eight modes used above (k y = 1,2,3, 
and 4 and a =| and [). 

We study eight samples of the system with different impurity configurations. As we mentioned, 
for each sample, the results depend sensitively on the initial state. Thus, we use the formula (3.8) 
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for each sample with the eight initial configurations (3.1). We obtain the conductance as a function 
of the width of the magnetic potential X w . In Fig.6, we compare the results of the present method 
with those obtained by the numerical treatment of the Kubo formula. The data are normalized 
by the conductance of the system without magnetic domain wall in order to see the effect of the 
domain wall. After the processing of (3.8), the data show roughly similar behavior. For example, 
if we look at the data for the sample of impurity- 1, the conductance decreases a little and then 
increases with the width. The sample of impurity-3 shows monotonic decrease with the width, etc.. 

The data in Fig. 6 depend strongly on the impurity configuration. We believe that this strong 
dependence is due to the smallness of the system. We find that the Kubo formula systematically 
gives larger normalized values, the reason why is under investigation now. 

§4. Quantum interference and conductance 

In this section, we study the effect of the width d x of the initial packet along the x direction. 
If the width is small, then the wave function is localized in the x direction and is spread out 
in momentum space. If the width is large then the wave packet is nearly a plane wave and the 
the quantum mechanical interference among spatially distant places is important. The width d x 
corresponds to the damping factor in the Kubo formula e (see Appendix A) in the sense that it 
causes decoherence in the system. 23 ) The small damping factor corresponds to small d x . 

In Fig. 7 and Fig. 8, we show the domain wall width dependence of the conductance for three 
values of e and d x , respectively. We find that the renormalized conductance is close to one in the 
strongly decoherent cases (Fig. 7(c) and Fig. 8(c)). This means that the effect of the domain wall 
becomes small as the decoherent effect increases. Thus we confirm that the intrinsic scattering by 
the magnetic potential requires a kind of quantum interference effect. Indeed, the increase of the 
conductivity in the presence of a magnetic domain wall has been interpreted as an effect of adiabatic 
motion (quantum mechanical coherent motion) of the electron in the magnetic potential. 8 ) 

§5. Effects of the magnetic domains on the conductance 

Making use of the present method, let us study the effect of magnetic domain wall on the 
conductance in a wire. The structure of the magnetization profile is given by 



Tatara and Fukuyama have calculated analytically the Kubo formula for the conductance by a 




(5.1) 




(5.2) 



perturbation theory. 8 ) They found 
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where L w = \/Dt w , D = % 2 k F T/2m 2 , r = Im/hkp, and t w = 24ir 2 \ w L x (p> b M ) 2 t/ (X f E f ) 2 . 

We study this problem by the wave packet method and a numerical evaluation of the Kubo 
formula. For this purpose we again use a periodic boundary condition and we set a pair of domain 
walls in the system. 

For the numerical calculation we use a 256 x 16 lattice, corresponding to a physical size of 
45Af x 2.8Af- Periodic boundary conditions are used. The impurity concentration is iii = 4.88%, 
V = 8.3E F , M = 0.2 ~ 0AE F , the temperature is set to E F /60, and e = 4.17 x 1(T 4 £ F . We 
study eight different impurity configurations and we plot the average of the results. The standard 
deviation of the data is about 0.1. The data strongly depend on the samples. This may be due 
to the smallness of the system. We find that the numerical data show the enhancement of the 
conductance, in agreement with the analytical results for the single domain wall, as shown in Fig.9. 

§6. Summary and Discussion 

We have studied the conductivity of systems with magnetic domain walls and impurity potential 
by direct numerical calculations of the time-dependent Schrddinger equation for the propagation 
of the electron wave packet (wave-packet method). We found that, when the system has magnetic 
domain walls, the transmission of the wave-packet depends sensitively on the initial configuration. 
We have proposed a method of a statistical treatment of data in the spirit of the multi-channel 
Landauer formula to obtain the electric conductivity from various initial configurations of the wave 
packets. In order to validate the method, we compare the results with the results obtained by direct 
numerical calculations of the conductivity by the Kubo formula. We conclude that the wave-packet 
method in combination with the statistical treatment is a useful approach to obtain information 
on the conductivity. 

We applied the method to a system with a ferromagnetic domain wall, and found that the 
analytical estimation of Kubo formula, numerical calculated Kubo formula and the wave-packet 
model consistently show the enhancement of conductivity in the magnetic domain wall. 

In the present work, we mainly concentrate on the comparison between the results of Kubo 
formula and the wave-packet method. A numerical evaluation of the Kubo formula requires a lot of 
computer memory because we use the simple L 3 method, where L is a size of the system. Memory 
and CPU time of the wave-packet method increase linearly with L and therefore it can be used for 
much large systems and for various boundary conditions. For instance, preliminary calculations for 
the three-dimensional case have been carried out. The results of the present study show sensitive 
dependence on the positions of the impurities. This may be attributed to the small size of the 
system whereas in large lattices we expect self-averaging of the conductivity. 

The authors would like to thank Dr. Gen Tatara for valuable discussions. The present work is 
partially supported by the Grant-in- Aid from the Ministry of Education. The calculation is done 
at the supercomputer center of ISSP. 



9 



Appendix: Direct calculation of the Kubo formula 

The dc conductivity can be obtained from the Kubo formula, 

CT0 = TT (£) / dE (-^)E(*|a')(a'|fc^(£-4)Mfi-^), (6.1) 

a,a' 

where \a) and E a are the a-th eigenstate and eigenvalue for the Hamiltonian, and f(E) is the 
Fermi distribution. 

In the Kubo formula approach, a non-zero current can be realized by adopting periodic boundary 
conditions. Eigenstates and eigenvalues are obtained by a direct numerical diagonalization of the 
system Hamiltonian. 

In a numerical approach, we must use a finite system, taking finite grids in real space, which leads 
to discretized energy levels. In this case, it is convenient to introduce the Lorenzian smoothing for 
the delta function in (6.1): 

5(x) ~ -ji-j, (6.2) 

with small e. In case of e — > 0, the conductivity vanishes due to the discreteness of the eigenvalues. 
In the present study, we take He to be the same order of energy spacing AE: 

he ~ AE. (6.3) 
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Fig. 1. Simulation model 
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Fig. 3. Probability distributions of the initial configurations for k y = 1, 2, 3 and 4 in Eq.(3.1) 
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Fig. 4. Dependence on the different initial configuration of the transmission coefficient as a function of the domain- 
wall width. Each symbol represents the data for different initial condition. 
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Fig. 5. Comparison of the transmission coefficients and conductivity. Kubo formula (solid line), Kubo formula with 
the effect of localization (dashed line), Kubo formula numerical method (closed circle), and the wave packet method 
with n — f,2,3, and 4 (cross, asterisk, open square, and closed square, respectively). 
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Fig. 6. Normalized conductance as obtained by the Kubo formula (a), and the wave-packet method (b). Different 
symbols represent data for different random potentials, using the same symbols in (a) and (b) for the same random 
potential. 
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Fig. 7. Dependence on e of the normalized conductance a/ao. (a):e = 4.17 x 10 _4 Ep,(b): e — 8.33 x lCT 4 EF,(c):e = 
4.71 x 10~ 3 Ep. Different symbols represent data for different random potentials, using the same symbols in (a), 
(b) and (c) for the same random potential. 
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Fig. 8. Dependence on d x of the normalized conductance a/ao- (a): d x = 3.54AF,(b): rf x = 2.12\F,(c):d x = 0.71Af. 
Different symbols represent data for different random potentials, using the same symbols in (a) , (b) and (c) for the 
same random potential. 
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Fig. 9. Comparison between the analytical (solid line) and the numerical results of the Kubo formula () and wave 
packet methodQ. 
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